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In this work, we find by means of first principle calculations a new physical mechanism to generate 
a two-dimensional electron gas, namely, the breaking of charge ordering at the surface of a charge 
ordered semiconductor due to the incomplete oxygen environment of the surface ions. The emergence 
of the 2D gas is independent of the presence of oxygen vacancies or polar discontinuities; this is 
a self-doping effect. This mechanism might apply to many charge ordered systems, in particular, 
we study the case of BaBiOa(OOl). In bulk, this material is a prototype of a forbidden valence 
compound in which the formal "metallic" Bi 4+ is skipped exhibiting a charge disproportionated 
Bi 3+ - Bi 5+ ordered structure. At room temperature, this charge disproportionation together with 
the breathing distortions gives rise to a Peierls semiconductor with monoclinic crystal structure. At 
higher temperatures (T> 750K) or upon doping, it turns cubic and metallic. Interestingly, doped 
BaBiOa was one of the first non-cuprates high-T c superconductors discovered. The outer layer of the 
Bi-terminated simulated surface turns more cubic- like and metallic while the inner layers remain 
in the insulating monoclinic state. On the other hand, the metallization does not occur for the Ba 
termination, a fact that makes this system appealing for nanostructuring. Finally, this finding sets 
another possible route for future exploration: the potential scenario of 2D superconductivity at the 
BaBiOa surface. 

PACS numbers: 73.20.-r, 71.45.Lr, 71.30.+h, 71.15.Mb 



The emergence of high mobility two-dimensional elec- 
tron gases (2DEGs) at the interfaces of insulating ox- 
ides has become a very exciting field of investigationjlj 
in the last decade. The reallization of 2DEGs at the 
oxide interfaces brings about a wide variety of phe- 
nomena (superconductivity [1 L magnetic-order Q, elec- 
tron correlation-driven effects 4j) that have awakened the 
interest on both fundamental issues and their future tech- 
nological applications. 

It is nowadays an issue of intense debate what is the 
origin of the 2DEG at the oxide interfaces. One point 
of view ascribes it to polar/non polar interfaces and is 
based on the polar "catastrophe" model that proposes 
an electronic reconstruction to compensate the growing 
dipole moment as the number of polar layers increases. 
Another invoked mechanism is the presence of oxygen va- 
cancies in the substrate. Each scenario explains part of 
the story and, probably, a complete understanding 
of the intrinsic nature of the 2DEG formation is highly 
dependent on the materials involved and the experimen- 
tal set up. Recently, it has been shown that a 2DEG can 
also be generated in a simpler context, namely, at the 
vacuum-cleaved surface of SrTi03. In this case, a metal- 
lic gas is formed independently of the oxide bulk carrier 
densities, opening the way towards novel means of 2DEG 
generation at the surface of transition-metal oxides Q. In 
this case, the presence of oxygen vacancies is suggested 
to lie behind the emergence of the metallic surface state. 

In this work, based on Density Functional 
Theory(DFT)-calculations0, we propose not only 



a new candidate able to sustain a confined electron gas 
but also a new physical mechanism to generate it that 
is different from the ones invoked until now. We show, 
namely, that a 2DEG is formed at the (001) surface of 
insulating BaBi03 as a consequence of a charge order 
disruption. No external factors, such as polar discon- 
tinuities or oxygen vacancies are necessary to obtain, 
in this case, the 2DEG except for the Bi- terminated 
surface itself. The surface generated carrier densities are 
quite high and of the same order of magnitude as the 
ones measured at the sharp LaAlOs/SrTiOa interfaces^ 
or at the cleaved SrTi03 surfaces. This phenomenon 
might be present in many other charge ordered materials 
as will be discussed later. 

In order to understand the nature of our finding we 
briefly describe the phase diagram of bulk BaBi03. At 
high temperature (T> 750 K) it is a cubic pcrovskite ex- 
hibiting metallic behaviour. Formally, one would expect 
each bismuth to have a valence 4+. However, Bi is a 
typical valence-skipping atom, and even in the high tem- 
perature metallic phase, it presents charge disproportion- 
ation. At lower temperatures, this disproportionation 
couples to the tilting of the BiOg octahedra. The crys- 
tal structure goes through a rhombohedral phase (750K 
> T > 405K), becoming monoclinic for T < 405K[uj. 
The charge disproportionation together with the struc- 
tural distortion are further enhanced giving rise to a 
formal Bi 5+ - Bi 3+ charge-ordered Peierls-likc insulator 
in the low temperature monoclinic phase. The oxy- 
gen octahedra around the Bi ions present alternating 
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FIG. 1: (Color online) Bandstructure and total density of 
states for: a) Ba- terminated (11 layers) and b) Bi- terminated 
(9 layers) BaBiOs (001). The corresponding DOS plots are 
inverted (energy vs DOS) and the units are in eV. 



brcathing-in and brcathing-out structural instabilities. 
Being the electronic properties of this material fascinant- 
ing by themselves, the discovery of high-T c superconduc- 
tivity in doped BaBiOaJTl], E3 made this system even 
more intriguing and interesting. It has been shown that 
upon doping, the monoclinic phase turns cubic or tetrag- 
onal (depending on the dopant) and metallic, exhibiting 
superconductivity with T c 's as high as 30K. 

Much of the understanding of the electronic structure 
and structural properties of both, the parent BaBi03 and 
the doped compounds, has been accomplished by means 
of first principle calculations [l3l - fl8| . In this material, the 
physics is dominated by spatially quite extended Bi(s)- 
O(p) orbitals. In the absence of important correlation ef- 
fects as in typical d or / electron systems, much progress 
has been made from early local-density approximation 
(LDA) calculations based on DFT. 

Calculations using LDA (or its gradient corrections 
GGA[2(|) could explain the splitting of the Bi(s)- 
O(p) band around due to Peierls- like distortions 
that are switched on, in particular, by the breath- 
ing instability [13]. These calculations could account, 
then, for the charge disproportionation and its relation 
to the structural distortions in the monoclinic phase. 
LDA(GGA) results predict a pseudomctallic behaviour, 



however, it is well known that BaBi03 bulk presents an 
indirect gap whose experimental reported values go from 
0.2 eV to 1.1 eV[2l|. More recently it has been shown 
that in order to open the indirect gap and describe quan- 
titatively the structural properties and the insulating be- 
haviour of this phase of BaBi03, it is necessary to go 
beyond standard DFT approaches, for instance by using 
hybrid functionals that combine a fraction of non-local 
exact exchange with local or semilocal aproximations [22^ . 

In this contribution, the theoretical study of the (001) 
surface of BaBiOa is faced for the first time. We per- 
form first principles DFT calculations and take care of 
the gap problem by cross-checking the results with func- 
tionals that go beyond LDA or GGA. All the results pre- 
sented in the main part of this work are done with a full- 
potential APW method in the GGA approximation as 
implemented in the Wien2k code[23j. In the Supplemen- 
tary section we show that the main physical findings are 
confirmed using more sophisticated functionals. The sur- 
faces are modeled by supercells with different slab thick- 
nesses. We consider slabs composed by 9 to 15 layers 
that are stacked following the monoclinic crystal struc- 
ture with both Ba- and Bi- terminated situations in the 
(001) direction. To avoid the interaction between oppos- 
site surfaces, they are separated in the z direction by an 
empty space volume ranging from 9 to 21 A. The super- 
cells have two inversion symmetric surfaces for simplicity. 
All internal atomic positions are allowed to relax. 

Our caculations indicate that Bi- terminated BaBiOs 
turns metallic while the Ba-terminated case remains in- 
sulating as in the bulk (monoclinic phase). The results 
obtained for different supercells with different number of 
layers are qualitatively the same in the corresponding ter- 
mination. The comparison among different slabs is useful 
to detect finite size effects (as described bellow). 

In Fig. [T]wc show the bandstructure and total densi- 
ties of states (DOS) for a Ba- and a Bi- terminated slab, 
a) and b), respectively. It can be clearly seen that in 
the first case, the system behaves as an insulator while 
in the second one there are several bands crossing E^, 
providing a considerably large amount of states visible 
in the corresponding DOS. These bandstructures corre- 
spond to slabs with 9 and 11 layers, having both of them 
ten Bi atoms in the slabs unit cell. In Fig. [T]a), the red 
bands have mainly Bi 3+ -0 character and the green ones 
mainly Bi 5+ -0 one. The blue bands have mostly O-p 
states with a significant Bi-6s weight. In Fig. [T]b), there 
are three green Bi 5+ -0 bands above E^ and three red 
Bi 3+ -0 bands below. One of the red bands crosses E^, 
presenting a small pocket at the M point. This crossing is 
due to a finite size effect. The contribution of this pocket 
to the Fermi surface decreases with growing slab thick- 
ness. There are finally, four black interpenetrating bands 
that drive the system metallic. They basically come from 
the surface Bi atoms which are strongly mixed with O-p 
states. 
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FIG. 2: (Color online) To the left, the simulated 13 layer slab with the Bi02 planes labeled from 1 to 4. The projected DOSs 
on these planes are in the central and right plots, for the Bi-s and O projected states, respectively. The botton row shows, for 
comparison, the corresponding DOSs for the bulk (monoclinic phase). Ef is at eV. 



The metallization of the Bi- terminated surface is a 
consequence of the incomplete octahedra environment 
of the Bi ions, which produces a rearrangement of the 
charge distribution suppresing the charge ordering at 
that BiC>2 plane. This suppression of the disproportion- 
ation is, in fact, partial but strong enough to turn the 
system metallic, as it is the case in the high temperature 
cubic phase of bulk BaBiC>3. On the other hand, for Ba 
surface termination, the oxygen octahedra environment 
of all Bi ions is complete, the charge ordering being then 
not affected in any of the BiC>2 planes. There is, indeed, 
a slight charge redistribution among the O atoms at the 
BaO surface plane, which has no effect on the insulating 
behavior of the whole system (24j|. 

In order to trace the origin of this metallicity, in Fig. [2] 
we plot the DOSs projected onto the Bi02 planes for the 
Bi- terminated 13 layer slab. The BaO layers are skipped 
for the sake of simplicity. The corresponding bulk pro- 
jected DOSs are presented in the bottom row of Fig. [2j 
for comparison. The effects of charge disproportionation 
in the bulk can be clearly observed within the [-1.3 eV, 
2 eV] energy range. There are quasi-symmetrically dis- 
tributed occupied and empty 6s bands around Hp, pre- 
senting the Bi 3+ ions mostly occupied states and the Bi 5+ 
ions mostly empty ones. As mentioned before, there is 
an important Bi(6s)-0(p) hybridization, which is appre- 
ciable in the O projected DOS. 

Our slab results show that already the third Bi layer 
from the surface, has a bulk like projected DOS. The 



slight downward shift of the Fermi level for layers 3 and 
4, as compared to the bulk, is again due to a finite size 
effect. The thinner the slab, the larger the downward 
shift of Ep into the 6s valence bands. The behaviour of 
the Bi-6s and O states of the surface and subsurface Bi- 
layers (labeled as 1 and 2) is qualitatively different to 
what happens in the deeper Bi planes. In layers 1 and 
2 there is an effective charge transfer from the originally 
Bi 3+ to the Bi 5+ ion. The system turns, in this way, 
metallic and this metallicity is mainly confined to the 
two outer Bi02 planes of the Bi- terminated BaBi03. 
For the thinner film considered, namely the 9 layer slab 
(not shown in Fig. [2]), the confinement is less effective 
but still the metallization is predominantly at the outer 
layers. 

The 2D carrier density, n2D, is obtained from Lut- 
tinger's theorem as n2D ~ vtti where Ap is the area 
enclosed by the 2D Fermi surface. The obtained n-io for 
the Bi-terminated slabs is around 0.6 electrons per 2D 
unit cell (or 1.6 x 10 14 cm -2 ). This carrier density is 
of the same order of magnitude as the ones estimated 
for sharp LaA103/SrTi03 interfaces and cleaved SrTi03 
surface 

The charge redistribution at the Bi-terminated sur- 
faces is present even in the unrelaxed systems, which turn 
metallic just by bond breaking. When the slabs are al- 
lowed to relax, the breathing distortions are washed out 
at the surface and subsurface Bi planes, contributing to 
an enhancement of the surface metallicity. In the deeper 
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FIG. 3: (Color online) Averaged Bi-0 bond lengths for the 
two types of Bi sites obtained for the 13 layer Bi-terminated 
BaBi03 slab. The label for each layer is the same as in Fig. 
El Squares (circles) correspond to the Bi 3+ -0(Bi 5+ -0) bond 
lengths, respectively. The solid and dashed straight lines are 
the corresponding values for low-T monoclinic and high-T 
cubic bulk phases, for comparison. 

layers, compressed and expanded octahedra remain with- 
out significant changes preserving there the insulating 
charge ordering as in the bulk. 

In Fig. [3] we plot the averaged Bi-0 bond length of 
the Bi 3+ and Bi 5+ ions of each layer for the 13 layer 
slab. We observe that the structure turns cubic like as 
one approaches the surface. There is experimental ev- 
idence from thin films measurements supporting these 
results (Hf. We can trace a parallelism between the effect 
of the surface on the structural and electronic properties 
of the Bi-terminated film and the effect of temperature on 
the same properties in bulk BaBi03. That is, we could 
think of having "cold" insulating monoclinic regions in 
the deeper planes and "hot" metallic cubic ones close to 
the surface. 

We can also make an analogy with the high-T c super- 
conductor, the doped BaBiC>3 in bulk, which turns cubic 
and metallic upon doping with K or Pb. In this context, 
the predicted 2DEG at the BiC>2 surface deserves further 
investigation regarding its superconducting properties. 

Another interesting finding with potential technolog- 
ical applications is the fact that the BaO surface is in- 
sulating. The possibility of drawing BiC<2 nanocircuits 
on top of BaO terminated BaBiC>3 surfaces constitutes a 
subject appealing for exploration. 

Finally, this surface metallization phenomenon might 
be present in many other charge ordered materiales. 
Potential candidates deserving further investigation are 
CaFe0 3 [26|, Pb 2 3 [27|, LuNi0 3 [H. 

Summarizing, we predict the formation of a 2DEG at 
the Bi- terminated surface of BaBiC>3 (001). The physi- 
cal mechanism behind the metallization of the surface is 
the breaking of charge order due to an incomplete oxygen 
octahedra environment of the Bi ions at the Bi02 surface. 
The obtained metallic state is confined to the outer lay- 
ers and presents a quite high 2D carrier density, of the 



order of 0.6 electrons per 2D unit cell. This phenomenon 
is probably not exclusive of BaBi03 and might ocurr in 
other charge ordered semiconductors. It is independent 
of any external factor such as the ambient oxygen pres- 
sure or polar discontinuities, making out of this system a 
self-doping surface with promising potential applications 
to oxide electronics. 
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I. SUPPLEMENTARY MATERIAL 



Local or semilocal functionals of the exchange-correlation potential, such as LDA or 
GGA, predict a pseudometallic behavior instead of the well known semiconducting nature 
of the low temperature phase of bulk BaBiCV In the main part of this contribution we 
predict, by means of DFT caculations using GGA, that a 2DEG is formed at the surface of 
a Bi-terminated BaBiO3(001) slab. We show there that the physical mechanism behind this 
metallization is a charge order breaking of the Bi ions at the surfuce. Although these GGA 



in 



calculations are not able to reproduce the insulating nature of the core of a sufficiently thick 
slab, they can very well account for the charge ordering disruption. It has been recently 
shown that to better describe the insulating behavior of bulk BaBiOa it is necessary to 
go beyong standard LDA or GGA functionals, for instance by using hybrid functionals that 
combine a fraction of non-local exact exchange with local or semilocal approximations. jH 2]. 

In the following we will show that the metallization of a Bi-terminated BaBi03 (001) 
slab is also obtained with more sophisticated functionals. We calculate the 9 layer slab 



using two other approachesr, the Heyd-Scuseria- 



and the modified Becke- Johnson potential (mB J) 
with the Vienna ab initio package (VASP) 



rnzerhof hybrid functional (HSE) 



5j. The HSE calculations are performed 



while the mBJ calculations are carried out 



with the Wien2k code [8j. The mBJ potential has been proposed by Tran and Blaha for a 
better description of band gaps. It is a local approximation to an atomic " exact-exchange" - 
potential and a screening term, plus the LDA correlation contribution. It has been shown 
that the mBJ potential provides a very good agreement with both the experimental data 
and the hybrid functionals or GW methods. The computational cost of the calculation 
using an mBJ potential is much less (by orders of magnitude) than the other two mentioned 
approaches, so that it can be applied to very large systems such as the ones we are interested 
here. 



As far as we know, the mBJ potential has not been tested on bulk BaBiOs yet. In Fig. [ST] 
we show the bandstructure and projected Bi-s DOS obtained by using (a) GGA, (b) mBJ, 
and (c) HSE. 

The opening of an indirect gap can be clearly seen in the results using the mBJ and 
HSE functionals as compared to the pseudometallic behavior obtained with GGA. The 
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FIG. SI: (Color online) Bandstructure and corresponding Bi-s projected DOS for bulk BaBiC>3 
using a) GGA, b) mBJ and c) HSE calculations. The DOS plots are inverted (energy vs DOS) and 
the units are in eV. The energy scale is similar to the one in Fig. 2. 

experimental reported values for this gap go from 0.2 eV to 1.1 eV j^j], so that both 
mBJ and HSE are well within this range. The position of the bands with mainly 2p 
character are also better re pro duced with mBJ and HSE when compared with experimental 
photoemission spectra jsl [lo| . Despite this gap problem, it can be noted in the corre- 
sponding projected Bi-s DOS, that the GGA results account reasonably well for Bi 3+ -Bi 5+ 
charge disproportionation. This is the key ingredient for the main results of this contribution. 



In Fig. IS2l we present the total DOS obtained with GGA, mBJ and HSE for the 9 layer 
Bi-terminated slab. We can see that the main 0-2p band moves to more negative energies 
following the sequence GGA-mBJ-HSE, similarly as in bulk BaBi03. Besides this difference, 
the metallization of the system in the three calculations is assertive. 

In Fig. [S3| we plot the bandstructure obtained using the mBJ potential. Around the 
Fermi energy, it looks very similar to the GGA one, presented in the Fig. lb) of the 
main part of the paper. Again, the main difference appears in the position of the 0-2p bands. 



The same layer by layer analysis of the projected DOS, as the one performed using GGA, 



2 




FIG. S2: (Color online) Total DOS for the 9 layer Bi- terminated BaBi0 3 (001) slab obtained with 
GGA, mBJ and 
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FIG. S3: (Color online) Bandstructure of the 9 layer Bi-terminated BaBi03 (001) slab. The color 
coding is the same as in Fig. 1 in the main part of the paper. 

applies for mBJ and HSE. The three calculations show that the metallic states are confined 
to the outer planes of the slab. 
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